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FIG. 91 
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FIG. 93 
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METHODS AND APPARATUSES FOR 
SIGNALING WITH GEOMETRIC 
CONSTELLATIONS IN A RALEIGH FADING 
CHANNEL 

5 

CROSS-REFERENCE TO RELATED 
APPLICATIONS 

The present invention claims priority as a Continuation-In- 
Part to U.S. patent application Ser. No. 12/156,989 entitled 10 
“Design Methodology and Method and Apparatus for Signal- 
ing with Capacity Optimized Constellation”, which issued as 
U.S. Pat. No. 7,978,777 and claims priority to U.S. Provi- 
sional Application Ser. No. 60/933,319 entitled “New Con- I5 
stellations for Communications Signaling: Design Method- 
ology and Method and Apparatus for the New Signaling 
Scheme” to Barsoum et al., filed Jun. 5, 2007. The present 
application also claims priority to U.S. Provisional Applica- 
tion Ser. No. 61/362,649 filed Jul. 8, 2010 entitled “Methods 20 
and Apparatuses for Signaling with Geometric Constellations 
in a Raleigh Fading Channel” to Barsoum et al. The disclo- 
sure of U.S. patent application Ser. No. 12/156,989 and U.S. 
Provisional Application No. 61/362,649 is expressly incor- 
porated by reference herein in its entirety. 25 

STATEMENT OF FEDERALLY SPONSORED 
RESEARCH 

This invention was made with Government support under 30 
contract NAS7-03001 awarded by NASA. The Government 
has certain rights in this invention. 

BACKGROUND 

35 

The present invention generally relates to bandwidth and/ 
or power efficient digital transmission systems and more spe- 
cifically to the use of unequally spaced constellations having 
increased capacity. 

The term “constellation” is used to describe the possible 40 
symbols that can be transmitted by a typical digital commu- 
nication system. A receiver attempts to detect the symbols 
that were transmitted by mapping a received signal to the 
constellation. The minimum distance (d„„„) between constel- 
lation points is indicative of the capacity of a constellation at 45 
high signal-to-noise ratios (SNRs). Therefore, constellations 
used in many communication systems are designed to maxi- 
mize d ml „. Increasing the dimensionality of a constellation 
allows larger minimum distance for constant constellation 
energy per dimension. Therefore, a number of multi-dimen- 50 
sional constellations with good minimum distance properties 
have been designed. 

Communication systems have a theoretical maximum 
capacity, which is known as the Shannon limit. Many com- 
munication systems attempt to use codes to be able to transmit 55 
at a rate that is closer to the capacity of a communication 
channel. Significant coding gains have been achieved using 
coding techniques such as turbo codes and LDPC codes. The 
coding gains achievable using any coding technique are lim- 
ited by the constellation of the communication system. The 60 
Shannon limit can be thought of as being based upon a theo- 
retical constellation known as a Gaussian distribution, which 
is an infinite constellation where symbols at the center of the 
constellation are transmitted more frequently than symbols at 
the edge of the constellation. Practical constellations are finite 65 
and transmit symbols with equal likelihoods, and therefore 
have capacities that are less than the Gaussian capacity. The 


2 

capacity of a constellation is thought to represent a limit on 
the gains that can be achieved using coding when using that 
constellation. 

Prior attempts have been made to develop unequally 
spaced constellations. For example, a system has been pro- 
posed that uses unequally spaced constellations that are opti- 
mized to minimize the error rate of an uncoded system. 
Another proposed system uses a constellation with equiprob- 
able but unequally spaced symbols in an attempt to mimic a 
Gaussian distribution. 

Other approaches increase the dimensionality of a constel- 
lation or select a new symbol to be transmitted taking into 
consideration previously transmitted symbols. However, 
these constellation were still designed based on a minimum 
distance criteria. 

SUMMARY OF THE INVENTION 

Systems and methods are described for constructing a 
modulation such that the constrained capacity between a 
transmitter and a receiver approaches the Gaussian chamiel 
capacity limit first described by Shannon [ref Shannon 1948]. 
Traditional communications systems employ modulations 
that leave a significant gap to Shannon Gaussian capacity. The 
modulations of the present invention reduce, and in some 
cases, nearly eliminate this gap. The invention does not 
require specially designed coding mechanisms that tend to 
transmit some points of a modulation more frequently than 
others but rather provides a method for locating points (in a 
one or multiple dimensional space) in order to maximize 
capacity between the input and output of a bit or symbol 
mapper and demapper respectively. Practical application of 
the method allows systems to transmit data at a given rate for 
less power or to transmit data at a higher rate for the same 
amount of power. 

One embodiment of the invention includes a transmitter 
configured to transmit signals to a receiver via a communica- 
tion channel, where the transmitter, includes a coder config- 
ured to receive user bits and output encoded bits at an 
expanded output encoded bit rate, a mapper configured to 
map encoded bits to symbols in a symbol constellation, a 
modulator configured to generate a signal for transmission 
via the communication channel using symbols generated by 
the mapper, where the receiver, includes a demodulator con- 
figured to demodulate the received signal via the communi- 
cation channel, a demapper configured to estimate likeli- 
hoods from the demodulated signal, and a decoder that is 
configured to estimate decoded bits from the likelihoods gen- 
erated by the demapper. In addition, the symbol constellation 
is a PAM-8 symbol constellation having constellation points 
within at least one of the ranges specified in FIGS. 25 - 47 . 

In a further embodiment, the code is a Turbo code. In 
another embodiment, the code is a LDPC code. 

In a still further embodiment, the constellation provides an 
increase in capacity over the Raleigh fading chamiel at a 
predetermined SNR that is at least 5% of the gain in capacity 
achieved by a constellation optimized for joint capacity at the 
predetermined SNR. 

In still another embodiment, the constellation provides an 
increase in capacity over the Raleigh fading chamiel at a 
predetermined SNR that is at least 20% of the gain in capacity 
achieved by a constellation optimized for joint capacity at the 
predetermined SNR. 

In a yet further embodiment, the constellation provides an 
increase in capacity at a predetermined SNR over the Raleigh 
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fading channel that is at least 50% of the gain in capacity 
achieved by a constellation optimized for joint capacity at the 
predetermined SNR. 

In yet another embodiment, the constellation provides an 
increase in capacity over the Raleigh fading channel at a 5 
predetermined SNR that is at least 90% of the gain in capacity 
achieved by a constellation optimized for joint capacity at the 
predetermined SNR. 

In another embodiment again, the constellation provides an 
increase in capacity over the Raleigh fading channel at a to 
predetermined SNR that is at least 1 00% of the gain in capac- 
ity achieved by a constellation optimized for joint capacity at 
the predetermined SNR. 

In a further additional embodiment, the constellation pro- 
vides an increase in capacity over the Raleigh fading channel 1 5 
at a predetermined SNR that is at least 5% of the gain in 
capacity achieved by a constellation optimized for PD capac- 
ity at the predetermined SNR. 

In another additional embodiment, the constellation pro- 
vides an increase in capacity over the Raleigh fading channel 20 
at a predetermined SNR that is at least 20% of the gain in 
capacity achieved by a constellation optimized for PD capac- 
ity at the predetermined SNR. 

In a still yet further embodiment, the constellation provides 
an increase in capacity over the Raleigh fading channel at a 25 
predetermined SNR that is at least 50% of the gain in capacity 
achieved by a constellation optimized for PD capacity at the 
predetermined SNR. 

In still yet another embodiment, the constellation provides 
an increase in capacity over the Raleigh fading channel at a 30 
predetermined SNR that is at least 90% of the gain in capacity 
achieved by a constellation optimized for PD capacity at the 
predetermined SNR. 

In still another embodiment again, the constellation pro- 
vides an increase in capacity over the Raleigh fading channel 35 
at a predetermined SNR that is at least 100% of the gain in 
capacity achieved by a constellation optimized for PD capac- 
ity at the predetermined SNR. 

A still further additional embodiment includes a transmit- 
ter configured to transmit signals to a receiver via a commu- 40 
nication channel, where the transmitter, includes a coder con- 
figured to receive user bits and output encoded bits at an 
expanded output encoded bit rate, a mapper configured to 
map encoded bits to symbols in a symbol constellation, a 
modulator configured to generate a signal for transmission 45 
via the communication channel using symbols generated by 
the mapper, where the receiver, includes a demodulator con- 
figured to demodulate the received signal via the communi- 
cation channel, a demapper configured to estimate likeli- 
hoods from the demodulated signal, and a decoder that is 50 
configured to estimate decoded bits from the likelihoods gen- 
erated by the demapper. In addition, the symbol constellation 
is a PAM-1 6 symbol constellation having constellation points 
within at least one of the ranges specified in FIGS. 48-71. 

Still another additional embodiment includes a transmitter 55 
configured to transmit signals to a receiver via a communica- 
tion channel, where the transmitter, includes a coder config- 
ured to receive user bits and output encoded bits at an 
expanded output encoded bit rate, a mapper configured to 
map encoded bits to symbols in a symbol constellation, a 60 
modulator configured to generate a signal for transmission 
via the communication channel using symbols generated by 
the mapper, where the receiver, includes a demodulator con- 
figured to demodulate the received signal via the communi- 
cation channel, a demapper configured to estimate likeli- 65 
hoods from the demodulated signal, and a decoder that is 
configured to estimate decoded bits from the likelihoods gen- 
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erated by the demapper. In addition, the symbol constellation 
is a PAM-32 symbol constellation having constellation points 
within at least one of the ranges specified in FIGS. 72-95. 

Another further embodiment includes a transmitter config- 
ured to transmit signals to a receiver via a communication 
channel, where the transmitter, includes a coder configured to 
receive user bits and output encoded bits at an expanded 
output encoded bit rate, a mapper configured to map encoded 
bits to symbols in a symbol constellation, a modulator con- 
figured to generate a signal for transmission via the commu- 
nication channel using symbols generated by the mapper, 
where the receiver, includes a demodulator configured to 
demodulate the received signal via the communication chan- 
nel, a demapper configured to estimate likelihoods from the 
demodulated signal, and a decoder that is configured to esti- 
mate decoded bits from the likelihoods generated by the 
demapper. In addition, the symbol constellation is a N-Di- 
mensional symbol constellation, where the constellation 
points in at least one dimension are within at least one of the 
ranges specified in FIGS. 25-95. 

BRIEF DESCRIPTION OF DRAWINGS 

FIG. 1 is a conceptual illustration of a communication 
system in accordance with an embodiment of the invention. 

FIG. 2 is a conceptual illustration of a transmitter in accor- 
dance with an embodiment of the invention. 

FIG. 3 is a conceptual illustration of a receiver in accor- 
dance with an embodiment of the invention. 

FIG. 4o is a conceptual illustration of the joint capacity of 
a channel. 

FIG. 4b is a conceptual illustration of the parallel decoding 
capacity of a channel. 

FIG. 5 is a flow chart showing a process for obtaining a 
constellation optimized for capacity for use in a communica- 
tion system having a fixed code rate and modulation scheme 
in accordance with an embodiment of the invention. 

FIG. 6 a is a chart showing a comparison of Gaussian 
capacity and PD capacity over the AWGN channel for tradi- 
tional PAM-2,4,8,16,32. 

FIG. 6b is a chart showing a comparison between Gaussian 
capacity and joint capacity over the AWGN channel for tra- 
ditional PAM-2,4,8,16,32. 

FIG. 7 is a chart showing the SNR gap to Gaussian capacity 
for the PD capacity and joint capacity over the AWGN chan- 
nel of traditional PAM-2,4,8,16,32 constellations. 

FIG. 8a is a chart comparing the SNR gap to Gaussian 
capacity of the PD capacity for traditional and optimized 
PAM-2,4,8, 16,32 constellations, over the AWGN channel. 

FIG. 8b is a chart comparing the SNR gap to Gaussian 
capacity of the joint capacity for traditional and optimized 
PAM-2,4,8, 16,32 constellations, over the AWGN channel. 

FIG. 9 is a chart showing Frame Error Rate performance of 
traditional and PD capacity optimized PAM-32 constellations 
in simulations involving several different length LDPC 
codes, over the AW GN channel. 

FIGS. 10(7-10(7 are locus plots showing the location of 
constellation points of a PAM-4 constellation optimized for 
PD capacity and joint capacity over the AWGN channel ver- 
sus user bit rate per dimension and versus SNR. 

FIGS. 11a and lib are design tables of PD capacity and 
joint capacity optimized PAM-4 constellations forthe AWGN 
channel in accordance with embodiments of the invention. 

FIGS. 12a-12rf are locus plots showing the location of 
constellation points of a PAM-8 constellation optimized for 
PD capacity and joint capacity over the AWGN channel ver- 
sus user bit rate per dimension and versus SNR. 
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FIGS. 13a and 13 b are design tables of PD capacity and 
joint capacity optimized PAM-8 constellations fortlie AWGN 
channel in accordance with embodiments of the invention. 

FIGS. 14a-14d are locus plots showing the location of 
constellation points of a PAM- 1 6 constellation optimized for 5 
PD capacity and joint capacity over the AWGN channel ver- 
sus user bit rate per dimension and versus SNR. 

FIGS. 15a and 15 b are design tables of PD capacity and 
joint capacity optimized PAM-16 constellations for the 
AWGN channel in accordance with embodiments of the to 
invention. 

FIGS. 16a-16d are locus plots showing the location of 
constellation points of a PAM-32 constellation optimized for 
PD capacity and joint capacity over the AWGN channel ver- 
sus user bit rate per dimension and versus SNR. 15 

FIGS. 17a and lib are design tables of PD capacity and 
joint capacity optimized PAM-32 constellations for the 
AWGN channel in accordance with embodiments of the 
invention. 

FIG. 18 is a chart showing the SNR gap to Gaussian capac- 20 
ity for traditional and capacity optimized PSK constellations 
over the AWGN channel. 

FIG. 19 is a chart showing the location of constellation 
points of PD capacity optimized PSK-32 constellations over 
the AWGN channel. 25 

FIG. 20 is a series of PSK-32 constellations optimized for 
PD capacity over the AWGN channel at different SNRs in 
accordance with embodiments of the invention. 

FIG. 21 illustrates a Q AM-64 constructed from orthogonal 
Cartesian product of two PD optimized PAM-8 constellations 30 
in accordance with an embodiment of the invention. 

FIGS. 22a and 22 b are locus plots showing the location of 
constellation points of a PAM-4 constellation optimized for 
PD capacity over a Raleigh fading channel versus user bit rate 
per dimension and versus SNR. 35 

FIGS. 23a and 23 b are locus plots showing the location of 
constellation points of a PAM-8 constellation optimized for 
PD capacity over a Raleigh fading channel versus user bit rate 
per dimension and versus SNR. 

FIGS. 24a and 24 b are locus plots showing the location of 40 
constellation points of a PAM- 16 constellation optimized for 
PD capacity over a Raleigh fading channel versus user bit rate 
per dimension and versus SNR. 

FIGS. 25-27 are tables showing the performance of geo- 
metric PAM-8 constellations optimized for Joint Capacity 45 
over a Raleigh fading channel at specific SNRs in accordance 
with embodiments of the invention. 

FIGS. 28-30 are tables listing the constellation points cor- 
responding to the geometric PAM-8 constellation designs 
optimized for Joint Capacity over a Raleigh fading chamiel at 50 
specific SNRs listed in FIGS. 25-27. 

FIGS. 31-32 are tables showing maximum ranges for the 
geometric PAM-8 constellation designs optimized for Joint 
Capacity over a Raleigh fading channel at specific SNRs 
listed in FIGS. 25-27. 55 

FIGS. 33-37 are tables showing the performance of geo- 
metric PAM-8 constellations optimized for PD Capacity over 
a Raleigh fading chamiel at specific SNRs in accordance with 
embodiments of the invention. 

FIGS. 38-43 are tables listing the constellation points cor- 60 
responding to the geometric PAM-8 constellation designs 
optimized for PD Capacity over a Raleigh fading channel at 
specific SNRs listed in FIGS. 33-37. 

FIGS. 44-47 are tables showing maximum ranges for the 
geometric PAM-8 constellation designs optimized for PD 65 
Capacity over a Raleigh fading chamiel at specific SNRs 
listed in FIGS. 33-37. 
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FIGS. 48-49 are tables showing the performance of geo- 
metric PAM-16 constellations optimized for Joint Capacity 
over a Raleigh fading chamiel at specific SNRs in accordance 
with embodiments of the invention. 

FIGS. 50-53 are tables listing the constellation points cor- 
responding to the geometric PAM- 16 constellation designs 
optimized for Joint Capacity over a Raleigh fading chamiel at 
specific SNRs listed in FIGS. 48-49. 

FIGS. 54-55 are tables showing maximum ranges for the 
geometric PAM-16 constellation designs optimized for Joint 
Capacity over a Raleigh fading chamiel at specific SNRs 
listed in FIGS. 48-49. 

FIGS. 56-59 are tables showing the performance of geo- 
metric PAM- 16 constellations optimized for PD Capacity 
over a fa Raleigh fading chamiel at specific SNRs in accor- 
dance with embodiments of the invention. 

FIGS. 60-67 are tables listing the constellation points cor- 
responding to the geometric PAM- 16 constellation designs 
optimized for PD Capacity over a Raleigh fading channel at 
specific SNRs listed in FIGS. 56-59. 

FIGS. 68-71 are tables showing maximum ranges for the 
geometric PAM- 16 constellation designs optimized for PD 
Capacity over a Raleigh fading chamiel at specific SNRs 
listed in FIGS. 56-59. 

FIGS. 72-73 are tables showing the performance of geo- 
metric PAM-32 constellations optimized for Joint Capacity 
over a Raleigh fading channel at specific SNRs in accordance 
with embodiments of the invention. 

FIGS. 74-81 are tables listing the constellation points cor- 
responding to the geometric PAM-32 constellation designs 
optimized for Joint Capacity over a Raleigh fading channel at 
specific SNRs listed in FIGS. 72-73. 

FIGS. 82-83 are tables showing maximum ranges for the 
geometric PAM-32 constellation designs optimized for Joint 
Capacity over a Raleigh fading channel at specific SNRs 
listed in FIGS. 72-73. 

FIGS. 84-85 are tables showing the performance of geo- 
metric PAM-32 constellations optimized for PD Capacity 
over a Raleigh fading chamiel at specific SNRs in accordance 
with embodiments of the invention. 

FIGS. 86-93 are tables listing the constellation points cor- 
responding to the geometric PAM-32 constellation designs 
optimized for PD Capacity over a Raleigh fading channel at 
specific SNRs listed in FIGS. 84-85. 

FIGS. 94-95 are tables showing maximum ranges for the 
geometric PAM-32 constellation designs optimized for Joint 
Capacity over a Raleigh fading channel at specific SNRs 
listed in FIGS. 84-85. 

DETAILED DESCRIPTION OF THE INVENTION 

Turning now to the detailed description of the invention, 
communication systems in accordance with embodiments of 
the invention are described that use signal constellations, 
which have unequally spaced (i.e. ‘geometrically’ shaped) 
points. In many embodiments, the communication systems 
use specific geometric constellations that are capacity opti- 
mized at a specific SNR. In addition, ranges within which the 
constellation points of a capacity optimized constellation can 
be perturbed and are still likely to achieve a given percentage 
of the optimal capacity increase compared to a constellation 
that maximizes & mi „, are also described. Capacity measures 
that are used in the selection of the location of constellation 
points include, but are not limited to, parallel decode (PD) 
capacity and joint capacity. 

In many embodiments, the communication systems utilize 
capacity approaching codes including, but not limited to, 
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LDPC and Turbo codes. As is discussed further below, direct 
optimization of the constellation points of a communication 
system utilizing a capacity approaching channel code, can 
yield different constellations depending on the SNR for 
which they are optimized. Therefore, the same constellation 5 
is unlikely to achieve the same gains applied across all code 
rates; that is, the same constellation will not enable the best 
possible performance across all rates. In many instances, a 
constellation at one code rate can achieve gains that cannot be 
achieved at another code rate. Processes for selecting capac- to 
ity optimized constellations to achieve increased gains based 
upon a specific coding rate in accordance with embodiments 
of the invention are described below. Constellations points for 
geometric PAM-8, PAM- 16, and PAM-32 constellations that 
are optimized for joint capacity or PD capacity at specific 15 
SNRs are also provided. Additional geometric PAM-8, PAM- 
16, and PAM-32 constellations that are probabilistically 
likely to provide performance gains compared to constella- 
tions that maximize d min , which were identified by perturbing 
the constellation points of geometric PAM-8, PAM- 16, and 20 
PAM-32 constellations optimized for joint capacity or PD 
capacity, are also described. The constellations are described 
as being probabilistically likely to provide performance 
gains, because all possible constellations within the ranges 
have not been exhaustively searched. Within each disclosed 25 
range, a large number of constellations were selected at ran- 
dom, and it was verified that all the selected constellations 
provided a gain that exceeds the given percentage of the 
optimal capacity increase achieved by the optimized constel- 
lations relative to a constellation that maximizes d mjn (i.e. a 30 
PAM equally spaced constellation). In this way, ranges that 
are probabilistically likely to provide a performance gain that 
is at least a predetermined percentage of the optimal increase 
in capacity can be identified and a specific geometric constel- 
lation can be compared against the ranges as a guide to the 35 
increase in capacity that is likely to be achieved. In a number 
of embodiments, the communication systems can adapt the 
location of points in a constellation in response to channel 
conditions, changes in code rate and/or to change the target 
user data rate. 40 

Communication Systems 

A communication system in accordance with an embodi- 
ment of the invention is shown in FIG. 1 . The communication 
system 10 includes a source 12 that provides user bits to a 
transmitter 14. The transmitter transmits symbols over a 45 
channel to a receiver 16 using a predetermined modulation 
scheme. The receiver uses knowledge of the modulation 
scheme, to decode the signal received from the transmitter. 
The decoded bits are provided to a sink device that is con- 
nected to the receiver. 50 

A transmitter in accordance with an embodiment of the 
invention is shown in FIG. 2. The transmitter 14 includes a 
coder 20 that receives user bits from a source and encodes the 
bits in accordance with a predetermined coding scheme. In a 
number of embodiments, a capacity approaching code such as 55 
a turbo code or a LDPC code is used. In other embodiments, 
other coding schemes can be used to providing a coding gain 
within the communication system. A mapper 22 is connected 
to the coder. The mapper maps the bits output by the coder to 
a symbol within a geometrically distributed signal constella- 60 
tion stored within the mapper. The mapper provides the sym- 
bols to a modulator 24, which modulates the symbols for 
transmission via the channel. 

A receiver in accordance with an embodiment of the inven- 
tion is illustrated in FIG. 3. The receiver 16 includes a 65 
demodulator 30 that demodulates a signal received via the 
channel to obtain symbol or bit likelihoods. The demapper 
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uses knowledge of the geometrically shaped symbol constel- 
lation used by the transmitter to determine these likelihoods. 
The demapper 32 provides the likelihoods to a decoder 34 that 
decodes the encoded bit stream to provide a sequence of 
received bits to a sink. 

Geometrically Shaped Constellations 

Transmitters and receivers in accordance with embodi- 
ments of the invention utilize geometrically shaped symbol 
constellations. In several embodiments, a geometrically 
shaped symbol constellation is used that optimizes the capac- 
ity of the constellation. In many embodiments, geometrically 
shaped symbol constellations, which include constellation 
points within predetermined ranges of the constellation 
points of a capacity optimized constellation, and that provide 
improved capacity compared to constellations that maximize 
d mirl are used. Various geometrically shaped symbol constel- 
lations that can be used in accordance with embodiments of 
the invention, techniques for deriving geometrically shaped 
symbol constellations are described below. 

Selection of a Geometrically Shaped Constellation 

Selection of a geometrically shaped constellation for use in 
a communication system in accordance with an embodiment 
of the invention can depend upon a variety of factors includ- 
ing whether the code rate is fixed. In many embodiments, a 
geometrically shaped constellation is used to replace a con- 
ventional constellation (i.e. a constellation maximized for 
d„„„) in the mapper of transmitters and the demapper of 
receivers within a communication system. Upgrading a com- 
munication system involves selection of a constellation and in 
many instances the upgrade can be achieved via a simple 
firmware upgrade. In other embodiments, a geometrically 
shaped constellation is selected in conjunction with a code 
rate to meet specific performance requirements, which can for 
example include such factors as a specified bit rate, a maxi- 
mum transmit power. Processes for selecting a geometric 
constellation when upgrading existing communication sys- 
tems and when designing new communication systems are 
discussed further below. 

Upgrading Existing Communication Systems 

A geometrically shaped constellation that provides a 
capacity, which is greater than the capacity of a constellation 
maximized for d mjll , can be used in place of a conventional 
constellation in a communication system in accordance with 
embodiments of the invention. In many instances, the substi- 
tution of the geometrically shaped constellation can be 
achieved by a firmware or software upgrade of the transmit- 
ters and receivers within the communication system. Not all 
geometrically shaped constellations have greater capacity 
than that of a constellation maximized for d min . One approach 
to selecting a geometrically shaped constellation having a 
greater capacity than that of a constellation maximized for 
d min is to optimize the shape of the constellation with respect 
to a measure of the capacity of the constellation for a given 
SNR. Another approach is to select a constellation from a 
range that is probabilistically likely to yield a constellation 
having at least a predetermined percentage of the optimal 
capacity increase. Such an approach can prove useful in cir- 
cumstances, for example, where an optimized constellation is 
unable to be implemented. Capacity measures that can be 
used in the optimization process can include, but are not 
limited to, joint capacity or parallel decoding capacity. 

Joint Capacity and Parallel Decoding Capacity 

A constellation can be parameterized by the total number 
of constellation points, M, and the number of real dimensions, 
N rflm . In systems where there are no belief propagation itera- 
tions between the decoder and the constellation bit demapper, 
the constellation demapper can be thought of as part of the 
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channel. A diagram conceptually illustrating the portions of a 
communication system that can be considered part of the 
channel for the purpose of determining PD capacity is shown 
in FIG. 4 a . The portions of the communication system that are 
considered part of the channel are indicated by the ghost line 
40 . The capacity of the channel defined as such is the parallel 
decoding (PD) capacity, given by: 


/-i 

C PD = J]l(X i -,Y) 

;=o 

where X, is the ith bit of the 1-bits transmitted symbol, and Y 
is the received symbol, and I(A;B) denotes the mutual infor- 
mation between random variables A and B. 

Expressed another way, the PD capacity of a channel can be 
viewed in terms of the mutual information between the output 
bits of the encoder (such as an LDPC encoder) at the trans- 
mitter and the likelihoods computed by the demapper at the 
receiver. The PD capacity is influenced by both the placement 
of points within the constellation and by the labeling assign- 
ments. 

With belief propagation iterations between the demapper 
and the decoder, the demapper can no longer be viewed as part 
of the channel, and the joint capacity of the constellation 
becomes the tightest known bound on the system perfor- 
mance. A diagram conceptually illustrating the portions of a 
communication system that are considered part of the channel 
for the purpose of determining the joint capacity of a constel- 
lation is shown in FIG. 4b. The portions of the communica- 
tion system that are considered part of the channel are indi- 
cated by the ghost line 42 . The joint capacity of the channel is 
given by: 

Joint capacity is a description of the achievable capacity 
between the input of the mapper on the transmit side of the 
link and the output of the channel (including for example 
AWGN and Fading channels). Practical systems must often 
‘ demap ’ channel observations prior to decoding. In general, 
the step causes some loss of capacity. In fact it can be proven 
that C G ^C J01NT s:Cp D . That is, C JOINT upper bounds the 
capacity achievable by C PD . The methods of the present 
invention are motivated by considering the fact that practical 
limits to a given communication system capacity are limited 
by C JOINT and C PD . In several embodiments of the invention, 
geometrically shaped constellations are selected that maxi- 
mize these measures. 

Selecting a Constellation Flaving an Optimal Capacity 

Geometrically shaped constellations in accordance with 
embodiments of the invention can be designed to optimize 
capacity measures including, but not limited to PD capacity 
or joint capacity. A process for selecting the points, and poten- 
tially the labeling, of a geometrically shaped constellation for 
use in a communication system having a fixed code rate in 
accordance with an embodiment of the invention is shown in 
FIG. 5 . The process 50 commences with the selection ( 52 ) of 
an appropriate constellation size M and a desired capacity per 
dimension r|. In the illustrated embodiment, the process 
involves a check ( 52 ) to ensure that the constellation size can 
support the desired capacity. In the event that the constellation 
size could support the desired capacity, then the process itera- 
tively optimizes the M-ary constellation for the specified 
capacity. Optimizing a constellation for a specified capacity 
often involves an iterative process, because the optimal con- 
stellation depends upon the SNR at which the co mmu nication 
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system operates. The SNR for the optimal constellation to 
give a required capacity is not known a priori. Throughout the 
description of the present invention SNR is defined as the 
ratio of the average constellation energy per dimension to the 
5 average noise energy per dimension. In most cases the capac- 
ity can be set to equal the target user bit rate per symbol per 
dimension. In some cases adding some implementation mar- 
gin on top of the target user bit rate could result in a practical 
system that can provide the required user rate at a lower SNR. 
10 The margin is code dependent. The following procedure 
could be used to determine the target capacity that includes 
some margin on top of the user rate. F irst, the code (e.g. LDPC 
or Turbo ) can be simulated in conjunction with a conventional 
15 equally spaced constellation. Second, from the simulation 
results the actual SNR of operation at the required error rate 
can be found. Third, the capacity of the conventional constel- 
lation at that SNR can be computed. Finally, a geometrically 
shaped constellation can be optimized for that capacity. 

20 In the illustrated embodiment, the iterative optimization 
loop involves selecting an initial estimate of the SNR at which 
the system is likely to operate (i.e. SNR,,,). In several embodi- 
ments the initial estimate is the SNR required using a con- 
ventional constellation. In other embodiments, other tech- 
25 niques can be used for selecting the initial SNR. An M-ary 
constellation is then obtained by optimizing (56) the constel- 
lation to maximize a selected capacity measure at the initial 
SNR,,, estimate. Various techniques for obtaining an opti- 
mized constellation for a given SNR estimate are discussed 
,(1 below. 

The SNR at which the optimized M-ary constellation pro- 
vides the desired capacity per dimension □ (SNR OIrt ) is deter- 
mined ( 57 ). A determination ( 58 ) is made as to whether the 
35 SNR„„, and SNR,„ have converged. In the illustrated embodi- 
ment convergence is indicated by SNR ow , equaling SNR,.,,. In 
a number of embodiments, convergence can be determined 
based upon the difference between SNR our and SNR,„ being 
less than a predetermined threshold. When SNR ou , and SNR,„ 
40 have not converged, the process performs another iteration 
selecting SNR ou , as the new SNR,„ ( 55 ). When SNR oii , and 
SNR,,, have converged, the capacity measure of the constel- 
lation has been optimized. As is explained in more detail 
below, capacity optimized constellations at low SNRs are 
45 geometrically shaped constellations that can achieve signifi- 
cantly higher performance gains (measured as reduction in 
minimum required SNR) than constellations that maximize 

TIWI* 

The process illustrated in FIG. 5 can maximize PD capacity 
50 or joint capacity of an M-ary constellation for a given SNR. 
Although the process illustrated in FIG. 5 shows selecting an 
M-ary constellation optimized for capacity, a similar process 
could be used that terminates upon generation of an M-ary 
constellation where the SNR gap to Gaussian capacity at a 
55 given capacity is a predetermined margin lower than the SNR 
gap of a conventional constellation, for example 0.5 db. Alter- 
natively, other processes that identify M-ary constellations 
having capacity greater than the capacity of a conventional 
constellation can be used in accordance with embodiments of 
60 the invention. For example, the effect of perturbations on the 
constellation points of optimized constellations can be used 
to identify ranges in which predetermined performance 
improvements are probabilistically likely to be obtained. The 
ranges can then be used to select geometrically shaped con- 
65 stellations for use in a communication system. A geometri- 
cally shaped constellation in accordance with embodiments 
of the invention can achieve greater capacity than the capacity 
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of a constellation that maximizes d min without having the 
optimal capacity for the SNR range within which the com- 
munication system operates. 

We note that constellations designed to maximize joint 
capacity may also be particularly well suited to codes with 
symbols over GF(q), or with multi-stage decoding. Con- 
versely constellations optimized for PD capacity could be 
better suited to the more common case of codes with symbols 
over GF(2) 

Optimizing the Capacity of an M-Ary Constellation at a 
Given SNR 

Processes for obtaining a capacity optimized constellation 
often involve determining the optimum location for the points 
of an M-ary constellation at a given SNR. An optimization 
process, such as the optimization process 56 shown in FIG. 5, 
typically involves unconstrained or constrained non-linear 
optimization. Possible objective functions to be maximized 
are the Joint or PD capacity functions. These functions may 
be targeted to channels including but not limited to Additive 
White Gaussian Noise (AWGN) or Rayleigh fading channels. 
The optimization process gives the location of each constel- 
lation point identified by its symbol labeling. In the case 
where the objective is joint capacity, point bit labelings are 
irrelevant meaning that changing the bit labelings doesn’t 
change the joint capacity as long as the set of point locations 
remains unchanged. 

The optimization process typically finds the constellation 
that gives the largest PD capacity or joint capacity at a given 
SNR. The optimization process itself often involves an itera- 
tive numerical process that among other things considers 
several constellations and selects the constellation that gives 
the highest capacity at a given SNR. In other embodiments, 
the constellation that requires the least SNR to give a required 
PD capacity or joint capacity can also be found. This requires 
running the optimization process iteratively as shown in FIG. 
5. 

Optimization constraints on the constellation point loca- 
tions may include, but are not limited to, lower and upper 
bounds on point location, peak to average power of the result- 
ing constellation, and zero mean in the resulting constella- 
tion. It can be easily shown that a globally optimal constella- 
tion will have zero mean (no DC component). Explicit 
inclusion of a zero mean constraint helps the optimization 
routine to converge more rapidly. Except for cases where 
exhaustive search of all combinations of point locations and 
labelings is possible it will not necessarily always be the case 
that solutions are provably globally optimal. In cases where 
exhaustive search is possible, the solution provided by the 
non-linear optimizer is in fact globally optimal. 

The processes described above provide examples of the 
manner in which a geometrically shaped constellation having 
an increased capacity relative to a conventional capacity can 
be obtained for use in a communication system having a fixed 
code rate and modulation scheme. The actual gains achiev- 
able using constellations that are optimized for capacity com- 
pared to conventional constellations that maximize d min are 
considered below. 

Gains Achieved by Optimized Geometrically Spaced Con- 
stellations 

The ultimate theoretical capacity achievable by any com- 
munication method is thought to be the Gaussian capacity, C G 
which is defined as: 


c G = yl°g 2 (l +snr ) 

Where signal -to-noise (SNR) is the ratio of expected signal 
power to expected noise power. Hie gap that remains between 
the capacity of a constellation and C G can be considered a 
measure of the quality of a given constellation design. 

10 The gap in capacity between a conventional modulation 
scheme in combination with a theoretically optimal coder can 
be observed with reference to FIGS. 6 a and 6b. FIG. 6 a 
includes a chart 60 showing a comparison between Gaussian 
capacity and the PD capacity of conventional PAM-2, 4, 8, 1 6, 
15 and 32 constellations that maximize d„„„. Gaps 62 exist 
between the plot of Gaussian capacity and the PD capacity of 
the various PAM constellations. FIG. 6b includes a chart 64 
showing a comparison between Gaussian capacity and the 
joint capacity of conventional PAM-2, 4, 8, 16, and 32 con- 
20 stellations that maximize d min . Gaps 66 exist between the plot 
of Gaussian capacity and the joint capacity of the various 
PAM constellations. These gaps in capacity represent the 
extent to which conventional PAM constellations fall short of 
obtaining the ultimate theoretical capacity i.e. the Gaussian 
25 capacity. 

In order to gain a better view of the differences between the 
curves shown in FIGS. 6 a and 6b at points close to the Gaus- 
sian capacity, the SNR gap to Gaussian capacity for different 
values of capacity for each constellation are plotted in FIG. 7. 
30 It is interesting to note from the chart 70 in FIG. 7 that (unlike 
the joint capacity) at the same SNR, the PD capacity does not 
necessarily increase with the number of constellation points. 
As is discussed further below, this is not the case with PAM 
constellations optimized for PD capacity. 

35 FIGS. 8a and 8b summarize performance of constellations 
for PAM -4, 8, 1 6, and 32 optimized for PD capacity and joint 
capacity (it should be noted that BPSK is the optimal PAM-2 
constellation at all code rates). The constellations are opti- 
mized for PD capacity and joint capacity for different target 
40 user bits per dimension (i.e. code rates). The optimized con- 
stellations are different depending on the target user bits per 
dimension, and also depending on whether they have been 
designed to maximize the PD capacity or the joint capacity. 
All the PD optimized PAM constellations are labeled using a 
45 gray labeling which is not always the binary reflective gray 
labeling. It should be noted that not all gray labels achieve the 
maximum possible PD capacity even given the freedom to 
place the constellation points anywhere on the real line. FIG. 
8 a shows the SNR gap for each constellation optimized for 
50 PD capacity. FIG. 8b shows the SNR gap to Gaussian capac- 
ity for each constellation optimized for joint capacity. Again, 
it should be emphasized that each '+’ on the plot represents a 
different constellation. 

Referring to FIG. 8a, the coding gain achieved using a 
55 constellation optimized for PD capacity can be appreciated 
by comparing the SNR gap at a user bit rate per dimension of 
2.5 bits for PAM-32. A user bit rate per dimension of 2.5 bits 
for a system transmitting 5 bits per symbol constitutes a code 
rate of 1/2. At that code rate the constellation optimized for 
60 PD capacity provides an additional coding gain of approxi- 
mately 1 .5 dB when compared to the conventional PAM-32 
constellation. 

The SNR gains that can be achieved using constellations 
that are optimized for PD capacity can be verified through 
65 simulation. The results of a simulation conducted using a rate 
1/2 LDPC code in conjunction with a conventional PAM-32 
constellation and in conjunction with a PAM-32 constellation 
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optimized for PD capacity are illustrated in FIG. 9. A chart 90 
includes plots of Frame Error Rate performance of the differ- 
ent constellations with respect to SNR and using different 
length codes (i.e. k=4,096 and k=16,384). Irrespective of the 
code that is used, the constellation optimized for PD capacity 5 
achieves a gain of approximately 1.3 dB, which closely 
approaches the gain predicted from FIG. 8 a. 

Capacity Optimized PAM Constellations 

Using the processes outlined above, locus plots of PAM 
constellations optimized for capacity can be generated that to 
show the location of points within PAM constellations versus 
SNR. Locus plots of PAM-4, 8, 16, and 32 constellations 
optimized for PD capacity and joint capacity and correspond- 
ing design tables at various typical user bit rates per dimen- 
sion are illustrated in FIGS. 10a-17fe. The locus plots and 15 
design tables show PAM-4, 8, 16; and 32 constellation point 
locations and labelings from low to high SNR corresponding 
to a range of low to high spectral efficiency. 

In FIG. 10a, a locus plot 100 shows the location of the 
points of PAM-4 constellations optimized for joint capacity 20 
plotted against achieved capacity. A similar locus plot 105 
showing the location of the points of joint capacity optimized 
PAM-4 constellations plotted against SNR is included in FIG. 

10 b. In FIG. 10c, the location of points for PAM-4 optimized 
for PD capacity is plotted against achievable capacity and in 25 
FIG. 10c? the location of points for PAM-4 for PD capacity is 
plotted against SNR. At low SNRs, the PD capacity optimized 
PAM-4 constellations have only 2 unique points, while the 
joint capacity optimized constellations have 3. As SNR is 
increased, each optimization eventually provides 4 unique 30 
points. This phenomenon is explicitly described in FIG. 11a 
and FIG. life where vertical slices ofFIGS. 10afe and 10cc? are- 
captured in tables describing some PAM-4 constellations 
designs of interest. The SNR slices selected represent designs 
that achieve capacities={0.5, 0.75, 1.0, 1.25, 1.5} bits per 35 
symbol (bps). Given that PAM-4 can provide at most log,(4) 

=2 bps, these design points represent systems with informa- 
tion code rates R={ 1/4, 3/8, 1/2, 5/8, 3/4} respectively. 

FIGS. 12afe and lied present locus plots of PD capacity 
and joint capacity optimized PAM-8 constellation points ver- 40 
sus achievable capacity and SNR. FIGS. 13a and 13fe provide 
slices from these plots at SNRs corresponding to achievable 
capacities r|={0. 5, 1.0, 1.5, 2. 0,2. 5} bps. Each of these slices 
correspond to systems with code rate R=q bps/log 2 (8), result- 
ing in R={l/6, 1/3, 1/2, 2/3, 5/6}. As an example of the 45 
relative performance of the constellations in these tables, 
consider FIG. 13fe which shows a PD capacity optimized 
PAM-8 constellation optimized for SNR=9.00 dB, or 1 .5 bps. 

We next examine the plot provided in FIG. 8a and see that the 
gap of the optimized constellation to the ultimate, Gaussian, 50 
capacity (CG) is approximately 0.5 dB. At the same spectral 
efficiency, the gap of the traditional PAM-8 constellation is 
approximately 1.0 dB. The advantage of the optimized con- 
stellation is 0.5 dB for the same rate (in this case R=l/2). This 
gain can be obtained by only changing the mapper and 55 
demapper in the co mmu nication system and leaving all other 
blocks the same. 

Similar information is presented in FIGS. 14a-14d, and 
15a-15fe which provide loci plots and design tables for PAM- 
16 PD capacity and joint capacity optimized constellations. 60 
Likewise FIGS. 16a-16<7, 17a and 17fe provide loci plots and 
design tables for PAM-32 PD capacity and joint capacity 
optimized constellations. 

Capacity Optimized PSK Constellations Over the AWGN 
Channel 65 

Traditional phase shift keyed (PSK) constellations are 
already quite optimal. This can be seen in the chart 180 
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comparing the SNR gaps of tradition PSK with capacity 
optimized PSK constellations shown in FIG. 18 where the 
gap between PD capacity and Gaussian capacity is plotted for 
traditional PSK-4, 8, 16, and 32 and for PD capacity opti- 
mized PSK-4, 8, 16, and 32. 

The locus plot of PD optimized PSK-32 points across SNR 
is shown in FIG. 19, which actually characterizes all PSKs 
with spectral efficiency r|<5. This can be seen in FIG. 20. 
Note that at low SNR (0.4 dB) the optimal PSK-32 design is 
the same as traditional PSK-4, at SNR=8.4 dB optimal PSK- 
32 is the same as traditional PSK-8, at SNR=14.8 dB optimal 
PSK-32 is the same as traditional PSK-16, and finally at 
SNRs greater than 20.4 dB optimized PSK-32 is the same as 
traditional PSK-32. There are SNRs between these discrete 
points (for instance SNR=2 and 15 dB) for which optimized 
PSK-32 provides superior PD capacity when compared to 
traditional PSK constellations. 

We note now that the locus of points for PD optimized 
PSK-32 in FIG. 19 in conjunction with the gap to Gaussian 
capacity curve for optimized PSK-32 in FIG. 18 implies a 
potential design methodology. Specifically, the designer 
could achieve performance equivalent or better titan that 
enabled by traditional PSK-4, 8, and 16 by using only the 
optimized PSK-32 in conjunction with a single tuning param- 
eter that controlled where the constellation points should be 
selected from on the locus of FIG. 19. Such an approach 
would couple a highly rate adaptive channel code that could 
vary its rate, for instance, rate 4/ 5 to achieve and overall (code 
plus optimized PSK-32 modulation) spectral efficiency of 4 
bits per symbol, down to 1/5 to achieve an overall spectral 
efficiency of 1 bit per symbol. Such an adaptive modulation 
and coding system could essentially perform on the optimal 
continuum represented by the rightmost contour of FIG. 18. 
Adaptive Rate Design 

In the previous example spectrally adaptive use of PSK-32 
was described. Techniques similar to this can be applied for 
other capacity optimized constellations across the link 
between a transmitter and receiver. For instance, in the case 
where a system implements quality of service it is possible to 
instruct a transmitter to increase or decrease spectral effi- 
ciency on demand. In the context of the current invention a 
capacity optimized constellation designed precisely for the 
target spectral efficiency can be loaded into the transmit map- 
per in conjunction with a code rate selection that meets the 
end user rate goal. When such a modulation/code rate change 
occurred a message could propagated to the receiver so that 
the receiver, in anticipation of the change, could select a 
demapper/decoder configuration in order to match the new 
transmit-side configuration. 

Conversely, the receiver could implement a quality of per- 
formance based optimized constellation/code rate pair con- 
trol mechanism. Such an approach would include some form 
of receiver quality measure. This could be the receiver’s 
estimate of SNR or bit error rate. Take for example the case 
where bit error rate was above some acceptable threshold. In 
this case, via a backchamiel, the receiver could request that 
the transmitter lower the spectral efficiency of the link by 
swapping to an alternate capacity optimized constellation/ 
code rate pair in the coder and mapper modules and then 
signaling the receiver to swap in the complementary pairing 
in the demapper/decoder modules. 

Geometrically Shaped QAM Constellations 

Quadrature amplitude modulation (QAM) constellations 
can be constructed by orthogonalizing PAM constellations 
into QAM in phase and quadrature components. Constella- 
tions constructed in this way can be attractive in many appli- 
cations because they have low-complexity demappers. 
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In FIG. 21 we provide an example of a Quadrature Ampli- 
tude Modulation constellation constructed from a Pulse 
Amplitude Modulation constellation. The illustrated embodi- 
ment was constructed using a PAM-8 constellation optimized 
for PD capacity for the AWGN channel at user bit rate per 
dimension of 1 .5 bits (corresponds to an SNR of 9.0 dB) (see 
FIG. 13 b). The label-point pairs in this PAM-8 constellation 
are {(000, -1.72), (001, -0.81), (010, 1.72), (011,-0.62), 
(100, 0.62), (101, 0.02), (110, 0.81), (111, -0.02)}. Exami- 
nation of FIG. 21 shows that the QAM constellation construc- 
tion is achieved by replicating a complete set of PAM-8 points 
in the quadrature dimension for each of the 8 PAM-8 points in 
the in-phase dimension. Labeling is achieved by assigning the 
PAM-8 labels to the LSB range on the in-phase dimension 
and to the MSB range on the quadrature dimension. The 
resulting 8x8 outer product fonns a highly structured QAM- 
64 for which very low-complexity de-mappers can be con- 
structed. Due to the orthogonality of the in-phase and quadra- 
ture components the capacity characteristics of the resulting 
QAM-64 constellation are identical to that of the PAM-8 
constellation on a per-dimension basis. The same process can 
be applied for constellations optimized for the Raleigh fading 
channel. 

N-Dimensional Constellation Optimization 

Rather than designing constellations in 1-D (PAM for 
instance) and then extending to 2-D (QAM), it is possible to 
take direct advantage in the optimization step of the addi- 
tional degree of freedom presented by an extra spatial dimen- 
sion. In general it is possible to design N-dimensional con- 
stellations and associated labelings. The complexity of the 
optimization step grows exponentially in the number of 
dimensions as does the complexity of the resulting receiver 
de-mapper. Such constructions constitute embodiments of 
the invention and simply require more ‘run-time’ to produce. 
Capacity Optimized Constellations for Fading Channels 

Similar processes to those outlined above can be used to 
design capacity optimized constellations for fading channels 
in accordance with embodiments of the invention. The pro- 
cesses are essentially the same with the exception that the 
manner in which capacity is calculated is modified to account 
for the fading channel. A fading channel can be described 
using the following equation: 

Y=ait)-X+N 


where X is the transmitted signal, N is an additive white 
Gaussian noise signal and a(t) is the fading distribution, 
which is a function of time. 

In, the case of a fading channel, the instantaneous SNR at 
the receiver changes according to a fading distribution. The 
fading distribution is Rayleigh and has the property that the 
average SNR of the system remains the same as in the case of 
the AWGN channel, E[X 2 ]/E[N 2 ]. Therefore, the capacity of 
the fading channel can be computed by taking the expectation 
of AWGN capacity, at a given average SNR, over a fading 
distribution of a, such as the Raleigh fading distribution, that 
drives the distribution of the instantaneous SNR. 

Many fading channels follow a Rayleigh distribution. 
FIGS. 22a-247) are locus plots of PAM-4, 8, and 16 constel- 
lations that have been optimized for PD capacity on a Ray- 
leigh fading channel. Locus plots versus user bit rate per 
dimension and versus SNR are provided. Similar processes 
can be used to obtain capacity optimized constellations that 
are optimized using other capacity measures, such as joint 
capacity, and/or using different modulation schemes. 
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Geometric PAM-8, PAM- 16, and PAM-32 Constellations 
As described above, geometric constellations can be 
obtained that are optimized for j oint or PD capacity at specific 
SNRs. In addition, ranges can be specified for the constella- 
5 tion points of a geometric constellation that are probabilisti- 
cally likely to result in geometric constellations that provide 
at least a predetermined performance improvement relative to 
a constellation that maximizes d„„„. Turning now to FIGS. 
25-95, geometric PAM-8, PAM- 16, and PAM-32 constella- 
to tions optimized for joint and PD capacity over the Raleigh 
fading channel at specific SNRs are listed. The performances 
of the optimal constellations are compared to the perfor- 
mances of traditional constellations that maximize d min . 
Ranges for the constellation points are also defined at specific 
15 SNRs, where constellations having constellation points 
selected from within the ranges are probabilistically likely 
(with probability close to one) to result in at least a predeter- 
mined performance improvement at the specified SNR rela- 
tive to a traditional constellation that maximizes d„„„. 

20 The geometric constellations disclosed in FIGS. 25-95 are 
defined by points y(i) such that y(i)=k(x(i)+r(i))+c. Values for 
x(i) and bounds on r(i) are provided in FIGS. 25-95 for 
PAM-8, PAM-16, and PAM-32 optimized for joint and PD 
capacity under Rayleigh fading channel conditions. For 
25 PAM-8 0<i<7, PAM- 16 0<i<15, and for PAM-32 0<i<31. To 
achieve optimal power efficiency, c should be set to zero. In 
addition to optimized constellations, FIGS. 25-95 specify 
ranges for the points of a geometric constellation, where 
selecting the points of a constellation from within the ranges 
30 is probabilistically likely to provide a geometric constellation 
having at least a predetermined performance improvement 
relative to a constellation that maximizes d min . The ranges are 
expressed as a maximum value for the constellation range 
parameter, r(i), which specifies the amount by which the point 
35 x(i) in the constellation is perturbed relative to the location of 
the corresponding point in the optimal constellation. A com- 
munication system using a constellation formed from con- 
stellations points selected from within the ranges specified by 
the maximum value (i.e. -r max sr(i)sr max ) is probabilistically 
40 likely to achieve a predetermined performance improvement 
relative to a constellation that maximizes d min . The predeter- 
mined performance improvements associated with the ranges 
specified in FIGS. 25-95 are expressed in terms of a percent- 
age of the increase in capacity achieved by the optimized 
45 constellation relative to a constellation that maximizes d m!M . 
Constellations formed from constellation points selected 
from within the ranges are probabilistically likely to achieve 
an increase in capacity at least as great as the indicated per- 
centage. 

50 With regard to the specific tables shown in FIGS. 25-95, 
each table is one of three different types of table. A first set of 
tables shows the performance of specific geometric constel- 
lations optimized for joint capacity or PD capacity. These 
tables include 6 columns. The first column enumerates a 
55 design number. The second column provides the SNR at 
which the constellation was optimized for the design defined 
by the entry in the first column. The third column provides the 
capacity achieved by the optimized constellation (Opt. Cap) 
at the SNR given in the second column. The fourth column 
60 provides the capacity achieved (Std. Cap) by a traditional 
uniformly spaced constellation i.e. a PAM constellation that 
maximizes d min (with the same number of points as the opti- 
mized constellation and where binary reflective gray labeling 
is assumed) at the SNR given in the second column. The fifth 
65 column shows the gain in bits per transmission provided by 
the optimized constellation over a constellation that maxi- 
mizes d,„,„. The sixth column shows the percentage gain in 
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capacity provided by the optimized constellation over the 
capacity provided by the traditional uniformly spaced con- 
stellation. 

A second set of tables lists the constellation points of the 
designs indicated in the first set of tables . These tables contain 
9 columns. The first column enumerates a designnumber. The 
remaining 8 columns describe a constellation point x(i) enu- 
merated by label in the second row of the table. Labels are 
given in decimal number format. With PAM 8 as an example, 
a label of 01 1 is given as the decimal number 3. 

The third set of tables specifies maximum perturbation 
ranges for the capacity optimized constellations indicated in 
the first set of tables, where the maximum ranges correspond 
to a high probabilistic likelihood of at least a predetermined 
performance improvement relative to a constellation that 
maximizes d min . These tables contain 8 columns. The first 
enumerates a design number (corresponding to a design from 
one of the aforementioned tables). The second column pro- 
vides the SNR for the design defined by the entry in the first 
column. The remaining 5 columns describe parameter x max 
which is the maximum amount any point in the designed 
constellation may be perturbed (in either the positive or nega- 
tive direction) and still retain, with probability close to unity, 
at least the gain noted by each column header of the joint or 
PD capacity as a percentage of the gain provided by the 
corresponding optimized point design over a traditional con- 
stellation that maximizes A min (all at the given SNR). Each 
table has a last column showing that if 100% of the gain 
afforded by the optimized constellation is desired, then 
parameter r(i) must be equal to zero (no deviation from 
designed points described in the point specification tables). 
Labelling of Constellations Using Cyclically Rotated Binary 
Reflective Gray Labels 

In performing optimization with respect to PD capacity, a 
conjecture can be made that constraining the optimization 
process to the subsequently described class of labelings 
results in no or negligible loss in PD capacity (the maximum 
observed loss is 0.005 bits, but in many cases there is no loss 
at all). Use of this labeling constraint speeds the optimization 
process considerably. We note that joint capacity optimiza- 
tion is invariant to choice of labeling. Specifically, joint 
capacity depends only on point locations whereas PD capac- 
ity depends on point locations and respective labelings. 

The class of cyclically rotated binary reflective gray labels 
can be used. The following example, using constellations 
with cardinality 8, describes the class of cyclically rotated 
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binary reflective gray labels. Given for example the standard 
gray labeling scheme for PAM-8: 

000,001,011,010, 110, 111, 101, 100 

Application of a cyclic rotation, one step left, yields: 

001,011,010, 110, 111, 101, 100,000 

Application of a cyclic rotation, two steps left, yields: 

011,010, 110, 111, 101, 100,000,001 

For a constellation with cardinality 8, cyclic rotations of 0 
to 7 steps can be applied. It should be noted that within this 
class of labelings, some labelings perform better than others. 
It should also be noted that different rotations may yield 
labelings that are equivalent (through trivial column swap- 
ping and negation operations). In general, labelings can be 
expressed in different but equivalent forms through trivial 
operations such as column swapping and negation operations . 
For example the binary reflective gray labels with one step 
rotation: 

001,011,010, 110, 111, 101, 100,000 

Can be shown to be equivalent to: 

000 , 001 , 011 , 111 , 101 , 100 , 110 , 010 

The above equivalence can be shown by the following steps 
of trivial operations: 

1 ) Negate the third column. This gives 

000 , 010 , 011 , 111 , 110 , 100 , 101,001 

2) Swap the second and third columns. This gives 

000,001,011,111,101,100,110,010 

The two labelings are considered equivalent because they 
yield the same PD Capacity as long as the constellation points 
locations are the same. 

In the constellation point specifications shown in FIGS. 
25-95, a labeling can be interchanged by any equivalent label- 
ing without affecting the performance parameters. A labeling 
used in the specifications may not directly appear to be a 
cyclically rotated binary reflective gray labeling, but it can be 
shown to be equivalent to one or more cyclically rotated 
binary reflective gray labelings. 

Prior .Art Geometric Constellations 
Geometric constellations have been specified in the prior 
art in attempts to achieve performance gains relative to con- 
stellations that maximize A min . Examples of such constella- 
tions are disclosed in Sommer and Fettweis, “Signal Shaping 
by Non-Uniform QAM for AWGN Channels and Applica- 
tions Using Turbo Coding” ITG Conference Source and 
Channel Coding , p. 81-86, 2000. The specific constellations 
disclosed by Sommer and Fettweis for PAM-8, PAM-16, and 
PAM-32 are as follows: 

PAM -8: 


- 1.6630 - 0.9617 - 0.5298 - 0.1705 0.1705 0.5298 0.9617 1.6630 

- 1.9382 - 1.3714 - 1.0509 - 0.8079 - 0.6026 - 0.4185 - 0.2468 - 0.0816 0.0816 


0.2468 0.4185 0.6026 0.8079 1.0509 1.3714 1.9382 


PAM-16: 


- 1.9382 - 1.3714 - 1.0509 - 0.8079 - 0.6026 - 0.4185 - 0.2468 - 0.0816 0.0816 

0.2468 0.4185 0.6026 0.8079 1.0509 1.3714 1.9382 


PAM-32: 


- 2.1970 - 1.7095 - 1.4462 - 1.2545 - 1.0991 - 0.9657 - 0.8471 - 0.7390 - 0.6386 

- 0.5441 - 0.4540 - 0.3673 - 0.2832 - 0.2010 - 0.1201 - 0.0400 0.0400 0.1201 

0.2010 0.2832 0.3673 0.4540 0.5441 0.6386 0.7390 0.8471 0.9657 

1.0991 1.2545 1.4462 1.7095 2.1970 
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Another class of geometric constellations is disclosed in 
Long et al., “Approaching the AWGN Channel Capacity 
without Active Shaping” Proceedings of International Sym- 
posium on Information Theory, p. 374, 1997. The specific 
PAM-8, PAM- 16, and PAM-32 constellations disclosed by 
Long et al. are as follows: 

PAM-8: 


-3 -1 -1 -1 1 1 13 


PAM-16: 
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constellation point (2.0691 , 1 .3438). The points correspond- 
ing to the remaining labels can be derived in a similar maimer. 

The ranges shown in FIGS. 25-95 can be utilized to select 
QAM constellations in a similar manner to that outlined 
5 above with respect to the selection of a PAM-8 constellation 
based upon ranges specified with respect to a PAM-8 constel- 
lation optimized for PD capacity at 9 dB. A range of 0.282 can 
be applied to every component of each two dimensional con- 
stellation point (based on the 50% of the achievable gain in 
10 capacity column). For example, the two points (-7.4992, 
-1.1438) and (2.2691, 1.1438) are within the ranges as they 
are spaced distances (-0.1, 0.2) and (0.2, -0.2) respectively 


-4 -2 -2 -2 -2 0 0 0 0 0 0 2 2 2 2 4 


PAM-32: 


-5 -3 -3 -3 -3 -3 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 

1 11 1 11 1 111333335 


The above prior art constellations are geometric and can 
provide performance improvements at some SNRs relative to 
constellations that maximize d mi „. The performance of the 
constellations varies with SNR and at certain SNRs the con- 
stellations are proximate to capacity optimized constella- 
tions. Therefore, the ranges specified in FIGS. 25-95 are 
defined so that prior art constellations are excluded at the 
specific SNRs at which these constellations are proximate to 
a capacity optimized constellation. 

Constructing Multidimensional Constellations 

The tables shown in FIGS. 25-95 can be used to identify 
optimal N-dimensional constellations. The optimized multi- 
dimensional constellation can be determined by finding the 
Cartesian power X" and the resulting labeling constructed by 
finding the corresponding Cartesian power of L". Ranges 
within which the multi-dimensional constellation points can 
be selected (i.e. perturbed), can then be defined with respect 
to each constellation point of the constructed multi-dimen- 
sional constellation, using an n-dimensional perturbation 
vector, such that each component of the perturbation vector 
has a magnitude that is less than r max defined by the range 
tables. 

Example of a QAM Constellation 

The optimized constellation points for a PAM-8 constella- 
tion optimized for PD capacity at SNR=9 dB are as follows: 


from the optimized constellation points. In this way, the 
ranges can be used to identify constellations that are proba- 
bilistically likely to result in a performance improvement 
relative to a constellation that maximizes d min . 

The same procedure can apply to a constellation optimized 
for joint capacity. However, the choice of labeling does not 
affect joint capacity. The above procedure can similarly be 
applied to an N-dimensional constellation constructed from a 
PAM constellation. 

Although the present invention has been described in cer- 
tain specific embodiments, many additional modifications 
and variations would be apparent to those skilled in the art. It 
is therefore to be understood that the present invention may be 
practiced otherwise than specifically described, without 
departing from the scope and spirit of the present invention. 
Thus, embodiments of the present invention should be con- 
sidered in all respects as illustrative and not restrictive. 

What is claimed is: 

1. A digital communication system, comprising: 
a transmitter configured to transmit signals to a receiver via 
a communication channel; 
wherein the transmitter, comprises: 

a coder configured to receive user bits and output 
encoded bits at an expanded output encoded bit rate; 
a mapper configured to map encoded bits to symbols in 
a symbol constellation; 


25 


30 


35 


40 


45 


- 7.3992 - 4.8128 - 1.3438 - 2.0696 7.3991 4.8129 1.3438 2.0691 


The labelings corresponding to the above PAM-8 constel- 
lation points are: 


000 001 010 011 100 101 no 111 


Using this PAM-8 constellation, it is possible to construct 60 
a QAM-64 constellation. While PAM-8 maps 3 bits to one 
dimension, QAM-64 maps 6 bits to two dimensions. The first 
three bits will determine the location in the X-dimension and 
the second three bits will determine the location in the Y-di- 
mension. The resulting QAM-64 constellation for example 65 
will map the bits 000 0 1 0 to the two dimensional constellation 
point (-7.3992, -1 .3438), and 111 1 10 to the two dimensional 


a modulator configured to generate a signal for transmis- 
55 sion via the communication channel using symbols 

generated by the mapper; 
wherein the receiver, comprises: 

a demodulator configured to demodulate the received 
signal via the communication channel; 
a demapper configured to estimate likelihoods from 
the demodulated signal; 

a decoder that is configured to estimate decoded bits 
from the likelihoods generated by the demapper; 
and 

wherein the points of the symbol, constellation have 
non-uniform inter-point distances selected to 
achieve a gain in capacity for a predetermined sig- 
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nal to noise ratio (SNR) relative to a symbol, con- 
stellation having uniform inter-point distances: and 
wherein the constellation provides an increase in 
capacity at a predetermined SNR that is at least 5% 
of the gain in capacity achieved by a constellation 
optimized for at least one of joint capacity and 
parallel decode PD capacity at the predetermined 
SNR in a Raleigh fading channel. 

2. The communication system of claim 1, wherein the code 
is a Turbo code. 

3. The communication system of claim 1, wherein the code 
is a LDPC code. 

4. The communication system of claim 1, wherein the 
constellation provides an increase in capacity at a predeter- 
mined SNR that is at least 20% of the gain in capacity 
achieved by a constellation optimized for joint capacity at the 
predetermined SNR in a Raleigh fading channel. 

5. The communication system of claim 1, wherein the 
constellation provides an increase in capacity at a predeter- 
mined SNR that is at least 50% of the gain in capacity 
achieved by a constellation optimized for joint capacity at the 
predetermined SNR in a Raleigh fading channel. 

6. The communication system of claim 1, wherein the 
constellation provides an increase in capacity at a predeter- 
mined SNR that is at least 90% of the gain in capacity 
achieved by a constellation optimized for joint capacity at the 
predetermined SNR in a Raleigh fading channel. 

7. The communication system of claim 1, wherein the 
constellation provides an increase in capacity at a predeter- 
mined SNR that is at least 100% of the gain in capacity 
achieved by a constellation optimized for joint capacity at the 
predetermined SNR in a Raleigh fading channel. 

8. The communication system of claim 1, wherein the 
constellation provides an increase in capacity at a predeter- 
mined SNR that is at least 20% of the gain in capacity 
achieved by a constellation optimized for PD capacity at the 
predetermined SNR in a Raleigh fading channel. 
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9. The co mmu nication system of claim 1, wherein the 
constellation provides an increase in capacity at a predeter- 
mined SNR that is at least 50% of the gain in capacity 
achieved by a constellation optimized for PD capacity at the 

5 predetermined SNR in a Raleigh fading channel. 

10. The communication system of claim 1, wherein the 
constellation provides an increase in capacity at a predeter- 
mined SNR that is at least 70% of the gain in capacity 
achieved by a constellation optimized for PD capacity at the 

10 predetermined SNR in a Raleigh fading channel. 

11. The communication system of claim 1, wherein the 
constellation provides an increase in capacity at a predeter- 
mined SNR that is at least 90% of the gain in capacity 
achieved by a constellation optimized for PD capacity at the 

15 predetermined SNR in a Raleigh fading channel. 

12. The communication system of claim 1, wherein the 
constellation provides an increase in capacity at a predeter- 
mined SNR that is at least 100% of the gain in capacity 
achieved by a constellation optimized for PD capacity at the 

20 predetermined SNR in a Raleigh fading channel. 

13. The communication system of claim 1, wherein the 
symbol constellation is a PAM-8 symbol constellation having 
constellation points within at least one of the ranges specified 
in FIGS. 25-47. 

25 14. The communication system of claim 1, wherein the 

symbol constellation is a PAM- 16 symbol constellation hav- 
ing constellation points within at least one of the ranges 
specified in FIGS. 48-71. 

15. The communication system of claim 1, wherein the 

30 symbol constellation is a PAM-32 symbol constellation hav- 
ing constellation points within at least one of the ranges 
specified in FIGS. 72-95. 

16. The communication system of claim 1, wherein the 
symbol constellation is a N-Dimensional symbol constella- 

3 5 tion, where the constellation points in at least one dimension 
are within at least one of the ranges specified in FIGS. 25-95. 



